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1 Introduction 


Our starting point in discussing gravitational lensing (hereafter GL) is Fermat’s po¬ 
tential i, i- Lens equation is given by the extremum of Fermat’s potential and the 
critical lines are points where the Jacobian from source plane to lens plane vanishes. 
(See Fig.l) 

Fermat’s potential is all from which we get information. Assuming thin lens ap¬ 
proximation, Fermat’s potential is given by Eq. (|2.1|) . Namely, Fermat’s potential (j) is 
hxed when we have determined the surface mass density Fl(a;) normalized by T^crit ■ 

(f) depends on the cosmological large scale structure through T,crit whereas Fl(a;) 
does on the lens models. The former problem is, for instance, concerned with the 
information of inhomogeneity of the universe, Dyer-Roeder distance, cosmological 
parameters, etc. The latter is related with the information of lensing object. 

Our arguments in this article are essentially restricted in the latter problem- 

model problem. 

As models, we may consider a variety of elliptical lens models that have different 
power dependence of distance in mass distribution. However their analytical sur¬ 
veys are not tractable and dedicated to numerical calculations. On the other hand, 
multipole expansion model seems to be rather tractable. 

In any case, the more models are sophisticated, the more parameter space gets 
complicated. Parameters are, for instance, mass, hnite core size, ellipticity, multipole 
moments and their angles, etc. Such parameters are hxed to reproduce the observed 
quantities such as image positions, image amplihcations and time delays between 
several images. However such parameter htting becomes ambiguous and best htting 
may not be unique in models with many parameters. Even if parameters had been 
htted uniquely, it does not mean systematic understanding of GL. 

More concretely speaking, observed data are biased in many cases by complicated 
unknown reasons. For instance, observed image amplihcations are disturbed by the 
dust lying between the source and observer. Also the magnitude of such disturbance 
may be diherent in the respective light path. In that case it is not necessarily reliable 
that we ht the parameters naively to the raw observed data. These problems are rather 
common to the observed cosmology. However they are especially serious in GL, since 
we have no conhrmed physical meaning in mass distribution model of lensing objects. 

So we need some strategy for searching the universal structure of GL that is 
irrelevant to models and boundary conditions. This is another problem — formulation 
problem. We apply catastrophe theory(CT) for such strategy. 

The application of CT to GL has already been performed in gravitational mi- 
crolensing and proved to be very useful |^. The concepts of optical depth, the time 
duration and amplihcation of image are made clear in the framework of GT. 

This paper is organized as follows. In §2 we argue general framework of singular¬ 
ities of caustics in complex coordinates, which is applied to the concrete models in 
subsequent sections. In §3 we discuss multipole expansion model. We will see how CT 
work in gravitational lensing. However it is indispensable to consider the transition of 
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singularities analytically for that purpose. These analyses are conhrmed by numerical 
calculations in §4 and applied to the real multiple quasar system. In the last section 
we consider elliptical lens model and see the consistency with the multipole expansion 
model. 
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2 General framework in complex coordinates 

Assuming the thin lens approximation, Fermat’s potential is given by 

(t){x^y) ='^{x — yY— — [ cPx'K{^)\n.\x — . (2.1) 

2 TT J 

Here x and y are two dimensional coordinates of image and source positions, 
respectively, which are illustrated in Fig.l . 


Fig.l 


K is the surface mass density of lens object. 


normalized by 




V = -Pq-s 

“ At^GDolDls ■ 

Thus we deal with the mapping of lens plane to source plane . 
So it becomes very convenient to introduce the complex coordinates. 


( 2 . 2 ) 


(ui, U2) = (u, u) = {xi + i X2 , Xi-i X2) ■ ( 2 . 3 ) 

Hereafter, we use Latin letters as the indices for real coordinates and Greek ones 
as the indices for complex coordinates. And we use the notation for derivatives, 
(pi = di(p = d(p/dxi, etc. Then we obtain the following expressions for the various 
quantities used to classify the singularities of caustics ||l|]. 


D = det{(Pij) = -A D = -Adet{<Pay) = A{\(Puu\^ - \(Puu\^ ) = [i - nf - |^|2 , (2.4) 

A(p = Tr(0ij) = A(puu = 2(1 - k) , (2.5) 

7 = 2(puu , (2.6) 

A = det(Zljj) = IGA = 16det(l)a/3) , (2.7) 

where 7 is the shear. Critical line is a curve in image plane characterized by H = 0 
where the amplihcation of images /r = T becomes 00 . Next the quantity that we 
often use in the classihcation of the singularities is the derivative operator along the 
critical line. 
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( 2 . 8 ) 


(f-V) = TA + ndr., 

where is a tangential vector to critical line, 

T, = . (2.9) 

So we obtain the expression 

L = (f ■ V) = ^{duDdu - duDdn) . (2.10) 

Here we can represent the conditions for various singularities in the complex co¬ 
ordinates. 

2.1 Singularities of Caustics 

1. Cuspoid sequence 

1-a. Fold singularity : 


O = 4.D = -[(l-fi;)2-|20„,|2] 
1-b. Cusp singularity : Besides 1-a condition, 

0 L(f)y^ — 7 {duD (l)uu duD (puu) 
2 i 

1-c. Swallowtail singularity ; Besides 1-b conditions, 

0 = 

1- d. Butterfly singularity : Besides 1-c conditions, 

0 = L^(j)u 

2. Umbilic : 

0a/3 = 0 {a,P = u,u) 

2- a. Elliptic Umbilic 

A > 0 

2-b. Hyperbolic Umbilic 


A < 0 


( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 
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2-c. Parabolic Umbilic 


3. Beak-to-beak and Lips 


A = 0 


3-a. Lips 


3-b. Beak-to-Beak 


duD = 0 


A > 0 


A < 0 


(2.17) 


(2.18) 


(2.19) 


2.2 Shapes of images 

Fermat’s potential is given by 


1 , 

= 2'“ 


u) , 


( 2 . 20 ) 


where is deflection potential and v = yi + i 1/2 represents sonrce position. Lens 
eqnation (Fermat’s principle) gives 


Differentiating 


V = u — 2dui^ = '^) ■ 

Eq-O), we obtain 


( 2 . 21 ) 


du ^ F 12 IF 12 ~ d^udv) ^ 1-[(1 - K)dv - -fdv] , 

I ^ u \ I I ^ 

= IF |2 ^ IF \2 ^-^ndv + F^dv) = ^[-^dv + { 1 - K)dv] , 

where 

= 1 - 2duduilJ = 1- K (real) , 

and 

F, ^ -2dl^ = 7 . 

Therefore Eq. (|2.4| ) is rewritten as 


( 2 . 22 ) 


(2.23) 
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So Icitip becomes 


\du\'^ = 


D = |FJ2 - |i^|2 . 


(1 — K)‘^\dv\‘^ r. ,, ,0 ,, dv 


D 2 


■[1 + |A| - (A—+ C.C.)] 


dv 


(2.24) 


(2.25) 


Here we have introduced the notation 


A = Fu/Fu = 


7 


1 — K 


(2.26) 


From Eqs. (|2.22|) and (|2.25|) , we know that the dv which gives the major axis of an 
image satisfies, 


arg(A|) = ^ , 

Therefore argument of du along major axis is 


(2.27) 


arg(dn) = ^arg( 7 ) + ^arg(l - «) + | 


(2.28) 


1 71 

= + 2^(1 - ’ 


(mod. vr) , 

where 6 is step function. From Eq.(|2.g|), the vector tangential to critical line is 


T = — D- 
2 % 

~^\,^u'4^uu T ‘^'4^uuu'4^UU T (1 ^)^'u] 


(2.29) 


Finally, we obtain the relative correlation between the direction of image elonga¬ 
tion and critical curve. 


nil 1 TT ~ 

arg(—) = -arg( 7 ) + -e{K - 1) - arg(T)s) 


(2.30) 


The meaning of Eq. (p.30|) is explained in the concrete models in subsequent sec¬ 
tions. Here we show how it works in the most simple model, spherical model. 

In this case, ip = 'p{\u\‘^) and 


7 = -2dlij = -2m V" 
where tp' = ^ with uj = 


etc. , 


(2.31) 


\U\ 


So we obtain the equation. 


1 TT 

axg{du) = axgu + -aigpip”) + 2 ~ ‘ 


(2.32) 
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If we assume that k is monotonically decreasing function of cn, we obtain the 
inequality, 


/ // 1 , K ,, 

(jJ L 

= [ dun] < 0 . (2.33) 

Jo 

Therefore the images are elongated along 

/ circular direction for k — 1 < 0 . 

\ radial direction for k — 1 > 0 . 


This conclusion is irrelevant to the image position relative to the critical lines. 



3 Multipole Expansion 


Multipole expansion is the expansion of Eq. (|2.1|) in powers of and is given by 


= vY ~ + 


m 


X ■ d 


X 


+ 


xqx 

1x1“^ 


(3.1) 


where m, d, q are rednced total mass, dipole and qnadrnpole moments of lens 
object, respectively. The ’’codimension” is identical with the nnmber of coordinates 
independent eqnalities. And the condition whether the singnlarity is generic or not is 
determined by Dimension of parameter space > Codimension. State variables x and 
control parameters y, m, d, q are expressed in complex coordinates as 


^ = xi + ix 2 , s = yi + iy 2 , 

d = di + id2 , q = qii + iqu + Q22 = 0 ) . 


(3.2) 


Taking, fnrthermore, the following rescaling. 


u , 

m2 


m2 


5 = 


d 


3 ) 
m2 


Q = 




we get Fermat’s potential in terms of complex coordinates 


(3.3) 


0 1, |2 11 I- I , , <5 

= — = 7 ^\u-v\ --In MM +-(- + -) + -(— + — 
m2 2 2 u u 2 


(3.4) 


The parameters in Eq.(p.4D are associated to physical parameters as follows: 
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M 

Merit 


Merit — 


Di + 0)2 

MeritDoi, 
III + ^12 

MeritDoL 


_ 3 

m 2 


m 


-2 






DolDos 

iGDLs 


(3.5) 

(3.6) 

(3.7) 


where A/, M ^rrel hx correspond the physical total mass, the physical 3-dimensional 
dipole moment and the physical 3-dimensional qnadrnpole moment, respectively. One 
of the aims of this article to hnd the relationships between model parameters and sin- 
gnlarities in analytical ways. So we hrst disenss abont the simplihed model and the 
lowest singnlarities. Onr task here is to express the eqnalities in §2 by the mnltipole 
parameters. 

Critical point condition Eq.( p.ll| ) becomes 

r, (5 Q 

m 2= l + 2- + 6^ . 3.8 

M M^ 


9 










Cusp condition Eg. ( p .121) is 

(1 + 2^ + 6^^^ 


\u\ 


(1 + 3 ^ + 12^)2 11 + 3 ^ + 12^|2 ■ 

V ?/. 11/^ / I 7/ 7/,^ I 


(3.9) 


Eqs.( p.8|) and ( |3.9| ) are too complicated for analytical surveys. So we discuss 6 and 
Q separately in the following. 

3.1 Q = 0 CASE 


At first, we consider Q = 0 case. Then Eqs. (|378|) and (|3.9|) are reduced to 

5, 


br =1 + 2- 


u 


\u\ 


From Eq. (|3.11|) , we obtain 


(i + zy _ 

(l + 3i)2 |l + 3ip' 


,(1 + 2 -)®, 

1ml . . j;, ] = 0 , 


(3.10) 

(3.11) 

(3.12) 

(3.13) 

where r is real and 6 is the physical angle of image position from the dipole 
direction. In terms of the geometrical variables, Eqs. (|3.101) , (|3.12|) and (|3.13|) are 
expressed by 


> 0 . 


■(l + 3^)2 

(1 + 2^ 

We introduce the geometrical variables. 




+ At cos 9 + A , 


sin6*[3t^ cos^ + t^(5 + 16 cos^ 9) + 48t cos6* + 36] = 0 , 


(3.15) 


(3.16) 


+ 12t‘^cos6* + 3t^(5 + 14cos^6*) -f 2t^ cos6*(51 + 16cos^6*) 

+12t(5 + 8 cos^ 0) + 72 cos 0 > 0 , 


(3.17) 


respectively, 
(i) sin 0 = 0 
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At first, we consider the critical points on the axis of the dipole direction. Then 
Eq. (p.l6|) vanishes. And Eqs. (p.lO|) and (|3.17|) are written as follows: 


|5|V = (l/ + 2r 


(3.18) 


l/(l/ + 2)>0, (3.19) 

where y = tcosO. The solutions of Eq. (|3.19|) change at |5| = a/3/9 qualitatively. 
At this value, beak-to-beak singularity appears. When |(5| is less than this value, cusp 
singularity appears. And when |5| is more than this value, fold singularity appears. 


(ii) sin 6* / 0 

Off the axis, we regard Eqs.( 3.1^ - (p.l7|) as the equations for S,tcos9{= y) and 
t^{= x). Using X and y, Eqs.( p.l5D and (|3.16| ) become 


| 5 | — X — 4: = Ay 


and 


X = 


4(2y + 3)- 
3?/+ 5 


respectively. However, x > 0,x > by dehnition, which implies 


Also Eq.( p.l7|) is rewritten as 


0 < 


-3<y < — 

_ y 3 


4x (y+ 2)^(y + 3) 
t (3?/+ 5)2 


Consequently y satishes 


-2 <y < — . 

^ 3 

Eliminating x from Eqs. (|3.20|) and (p.21|) , we obtain 

A\5\\2y + 3f = {3y + h)\y + 2\ 

That is. 


2 9 p(l - p) _ 




4(2p+ 1)3 


= /(P) 


where p is bounded 


(3.20) 

(3.21) 


(3.22) 

(3.23) 

(3.24) 

(3.25) 


0 < p = — (3t cos 6* + 5) < 1 . 
The function /(p) relation is plotted in Fig.2. 
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Fig. 2 


Two solutions coalesce at |(5p = in which higher singularity, swallowtail, 

appears. At |(5| more (less) than this value, cusp (fold) singularity appears. 


3.2 6 = 0 CASE 


Next, we consider 6 


0 case. Then Eqs. (|378D and (p^ are 


, ,9 , Q, 

Mp = 1 + 6-^ , 

(3.26) 

(l + 6g)3 _ |n|6 

(1 + 125)2 |1 + 125|2 • 

(3.27) 

From Eq.(|3.27) we obtain 


. r 1 n 

■“‘(1 + 12 +^’=°' 

(3.28) 

(1 + +)’ . „ 

(1 + 12 S )" 

(3.29) 

We introduce the geometrical variables 


7/2 t2 

^ ' 2?/9 

Q IQI 

(3.30) 


where r is real and 6 is the physical angle of image position from the quadrupole 
direction. In terms of the geometrical variables, Eqs. (p.26| ), (|3.28| ) and ( |3.29|) are 
expressed by 


\Q\w'^ = y/w'^ + 12w cos 29 + 36 , (3.31) 

sin 26[w‘^ + 12w^ cos 26 — 36ta^(4 cos^ 29 — 1) 

-8 ■ Q^w cos 20 - 4 ■ 6^] = 0 , (3.32) 


+ A2w^ cos 29 + 36ic'^(5 + 14 cos^ 29) + cos 20(21 + 4 cos^ 20) 

+8 ■ + cos^ 20) + 4 ■ 96w cos 20 > 0 , (3.33) 


respectively. 
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(i) sin 29 = 0 

At first, we consider the singularity on the axis of the quadrupole direction and 
on the axis normal to it. Then Eqs.(3.31) and (|3.33|) are written as follows: 


\Q\Y = iy + Gy 


y{y + 6)>0, (3.34) 

where y = |m/QPcos 26'. The solutions of Eqs. (|3.34|) change at |Q| = 1/24 quali¬ 
tatively. At this value, beak-to-beak singularity appears. When |Q| is less than this 
value, cusp singularity appears. And when |Q| is more than this value, fold singularity 
appears. 


(ii) sin 29 ^0 

we regard Eqs. ( |3.31|) - ( |3.33|) at points off those axes as the equations for Q, w cos 29 

and Those solutions of Eqs. (|3.31|) - ( |3.33|) qualitatively change at |Q| = 
and IQI = 1/8. They correspond to swallowtail and butterfly singularity, respectively. 
When IQI is less than the swallowtail-value, cusp singularity appears. And when IQI 
is more than that value, fold singularity appears. Talking about butterfly singularity, 
when IQI is not butterfly-value, cusp singularity always appears. These results are 
conhrmed by numerical calculations. 


3.3 Shapes of images 

From Eq.(p.9|) the vector tangential to critical line is 


Tu ^.'4^uu'^uuu 

2i 8i 

(3.35) 

for non-transparent lens model like multipole expansion. 
(2.30|) it follows 

From Eqs. (13.351) and 

arg(^) = ^arg(7) - argid^j) 


2 2 

(3.36) 

If source position is in the neighborhood of cusp singularity, is concluded 

to be almost negative dehnite from Eq.(|3.9|). Therefore, 


du 

arg(—) ~ 0 or tt . 

U 

(3.37) 


Namely, the images in the neighborhood of cusp are parallel to critical line. 
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For the general source position, 


arg( 


du, 

■L 1/ 


3,50, r 5 g, 

= T^arg 1 + 2- + 6^ - arg 1 + 3- + 12-^ 
2 u u 


(3.38) 
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4 Numerical Calculations 

In the previous section we have analyzed lens equation in which lensing object is 
modeled by multipole expansion. In this section we will check the results of each case 
by numerical calculations in the first two subsections. In the last subsection we will 
apply multipole expansion model to PG1115+080 §], which is compared with the 
article by Yoshida and Omote [Q . 

4.1 Q = 0 case 

Fig.3 shows the behaviors of caustics and critical lines in the neighborhood of |(5| = ^ 
at which beak-to-beak appears on axis. 


Fig. 3 


In |h| < the inequality in Eqs. (|3.19|) is satisfied and cusp appears on axis 
(Fig.3a). At |5| = ^ cusps coalesce to beak-to-beak on axis (Fig.3b). In |(5| > 
the inequality in Eqs.( |3.11^ ) is not satished. So cusps disappear and fold is born 
(Fig.3c). 


Fig.4 shows the behaviors in the neighborhood of |5p = 
appears off axis. 


2 _ 10+7vA 


63 


at which swallowtail 


Fig.4 


This corresponds to the argument in subsection §3.1(ii). As was indicated there, 

cusp singuralities at points P change to swallowtail when |(5| increases to 6c = \J 
(Fig.4a and Fig.4b) and cusps disappear and fold appears(Fig.4c). Swallowtail singu¬ 
larity is easily checked by Fig.4a corresponding to |5| < 5c- 


4.2 |(5| = 0 Case 

First we simulate the case (i) in §3-2. Fig.5 shows the behaviors of singularities in the 
neighborhood of |Q| = ^ at which beak-to-beak appears on 0 = | axis. 


Fig. 5 
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At IQ I < the inequality in Eqs. (|3.3^) is satisfied and cusps appear (Fig.5a). 
They coalesce to beak-to-beak at |Q| = A on the F-axis (Fig.5b). The inequality in 
Fqs. (|3.34|) is not satisfied at |Q| > A and fold appears instead of cusp (Fig.5c). 
Furthermore the case (ii) is simulated in Fig.6 and Fig. 7. 


Fig. 6 


Fig. 7 


IQI increases passing through two critical points, |(= Qi) and Q 2 ). 

Fig.6b and Fig.7b correspond to the transition points of butterfly and swallowtail, 
respectively. This is confirmed from Fig.6c and Fig.7a that are diagrams with |Q| > 
Qi and |Q| < Q 2 , respectively. 

4.3 Application to multiple quasar 

In this subsection we apply multipole expansion model Fq. ( |3.4| ) to lensing system 
PGl 115+080. We adopt the observed data by Christian et al 0, which is exhibited 
in Table 1. 


Table 1 


Numerical result is shown in Fig.8. 


Fig. 8 


Here we have taken the values for control parameters in Fq. (|3.4| ) as follows: 

source position = (O'kOOO, 0'7021) , 


m = 3.468 X 10“^^ , 


(4.1) 


1^1 = 0.1912 , 9s = 91°.85 , 

IQI = 0.0864 , 9q = 135°.79 , 
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where 6s and 6q are arguments of 6 and Q, respectively. The parameters adopted 
in are related to our parameters by 


= m (mass) , 

1^1 

= m^/^\6our\ = 1.126 X 10“® , 

lel 

= 2m|g| = 5.993 x lO'^^ , 

Xd 

0 

= es--j = 23°.95 , 

Xg 

Oq 

= AL = 67°.90 . 

2 


(4.2) 

Calculated image positions are 


e{Ai) = (-0".84, -0".87) , 
e{A2) = (-l".12,-0".26) , 
e{B) = (0".70, -0".60) , 
e(C) = (0".35,1".34) , 
e{D) = (0".56, -0".42) , 
e{E) = (-0".49,0".08) , 


(4.3) 


where we had the additional two images labeled by D and E as was indicated by 
Ref.[Q. Here we show source position and image positions 6 in unnormalized values 
that can be compared with observed values, while we use variables u, v normalized 
by ^/m in Eq. (^(4|) . 

The relative magnihcations of images are also calculated from Eq. (|2.4|) as 


1 

= 2.55 , 

1 6'A2 I 

= 2.22, 1^ 

fdC 

6-0 

6c 

/^D 

= 0.15 , 

]^6e 

= 0.73 X 10-2 


6c 



(4.4) 


We will check the validity of approximation in multipole expansion. It must hold 
the following conditions 


, 1^1 \Q\ 

1 > -^4 > ^ 


\u\ 


The best £t parameters give 

1^1 


1 —r 0.2266 and 

\u\ 


\u\ 


IQI 


\u\ 


0.1390 


(4.5) 


(4,6) 


Thus the condition Eq.(^) is rather marginal. So we should check that 2^-pole 
and higher multipole terms less contribute, though we do not touch this problem in 
this paper. 
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5 Discussions 


We have argued the classification of singularities of caustics based on the multipole 
expansion model of lensing object. In the application to the multiple quasar system of 
PG1115+080 there have arisen several problems. Fig .8 shows that deamplified images 
D and E are both normal to critical line. Arg(^) is given by 


from Eq. ( p.30| ) and k 
by 


arg(^) = ^arg( 7 ) 


arg(79s7) 


(5.1) 


= 0 for multipole expansion. Whereas amplification is given 




(5.2) 


i-i7r 

Amplification does not depend on the argument of 7 . 

Therefore if we are concerned with fold catastrophe we can not relate the ampli¬ 
fication with the image direction relative to the critical line. 

We should also be careful to the validity of approximation. Image D and E are 
near to the lens position and multipole expansion is badly convergent. We pointed 
out that dipole and quadrupole terms give the contribution of same order. 

It is certainly desirable to consider the same lensing system based on the several 
models and to check the consistency of models. 

We will compare our results with those of elliptical lens model by Narasimha et 
al p. Lensing object in p is modeled as the oblate spheroid 


+ y‘^ + 


l-e2 


= a 


(5.3) 


Mass density is given by 


p{a) = po(l + ^) ^ for ajrc<n 
ro 

= 0 for a/vc > n , (5.4) 

where Tc is core size. We can always set x,y coordinates in lens plane and z-axis 
is tilted relative to the photon path by an angle (p. 

Then the surface mass density projected on the lens plane is 


where 


m = 


1 — e 


2 2 _ ’ 


1 — sin (p 


(5.5) 






(5.6) 
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and 


6 ^ = + 


y 


(5.7) 


1 — sin^ tp 

In this model the dipole component, 5, vanishes and we will estimate the quadrupole 
component, adopted in P as 


e = 


1 / dxT,{x)z^ 

Dl^ /dfS(f) 
sin^ p db b^Kk{b) 
2Dl^ J^dbbK.ib) 


sin^ p rl 


where z is given by Eq. 


2 Dl^ Ff) ’ 

and where is defined by 


(5.8) 


Fk^ = f dx [ dy ^ - (1 + 1/) 

JO Jx Jy — X 


Vy 


[ dxx''{l + x) 2 [ du{l + u^) ■■ 
Jo Jo 


As is easily checked 


and 


Fg^^ = 2[ln(A/rF~TT + n) 


n 


a/?F”TT 


(5.9) 


(5.10) 


Fi°) + Fi^^ = - 


2 


Therefore we obtain 


F 


( 1 ) 


3 + 1 


48.5 for n=20 


( 6 . 11 ) 


(5.12) 


p(0) “ 95.9 for n=30 . 

-^3 t 

rdD ql in Eq.( p.8|) is related to the numerical result by 

^ = h£5l), (5.13) 

Dql S 

where ^(rad) is the observed image separation angle in unit of radian and s is the 
separation length in unit of core size. 

So 
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(5.14) 


2 ' s 

We adopt the values from Table 1 in 


e^sin^ip ^(sec) ^ 4 in-ii 

^ —-) ^77^ X 2.4 X 10 


40 ) 


AGMD 

k = 3 , esino9 = 0.6 , u = 20 , ——— = 32 , and sbc = 7.46 , (5.15) 




where sbc s between images B and C. Combining these values with 6bc — 1^7989 
we obtain 

^ = 1.51 X 10“^^ . (5.16) 

If we adopt other values ^ takes the range 


^elliptical — (2 ~ 4) X ^multipole 


(5.17) 


This seems that quadrupole moments evaluated by the different models coincide 
in this order. 

We may set fc = 1 or fc = 5 in elliptical lens model. Especially k = 2 seems 
to be important for isothermal model 0|]T3- In this case, numerical calculation 


becomes cumbersome and will be discussed elsewhere. Lens models so far considered 
are somehow empirical. By applying the several models into the same lensing system, 
we can construct the more reliable model of lensing object and gravitational lensing. 

We could not still exhaust the merit of CT in this article. In the near future, the 
numbers of observed lensing events will rapidly increase and we will get more precise 
information [O. CT will play increasingly essential role in that situation. 
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Figure Captions 


Fig.l Schematic diagram of the gravitational lens geometry. The solid line shows the 
real path of the light ray. 

Fig.2 Graph of the function /(p) given in Eq.( p.25|) . Two solutions of /(p) = 
coalesce at |5p = ^ maximum of /(p). 

Fig.3 Caustic and critical lines in neighborhood of |5| = a/3/9 at which beak-to-beak 
appears on axis. 

Fig.3a |(5| = 0.17 
Fig.3b 1^1 = 73/9 ~ 0.19 
Fig.3c |(5| = 0.21 

Fig.4 Caustic and critical lines in neighborhood of |5| = 7 = which 

swallowtail appears. Points P in Fig.4a are the points of cusp singularities 
which change to swallowtail when |(5| increases. 

Fig.4a |(5| = 0.32 
Fig.4b |5| = 5c ~ 0.36 
Fig.4c |5| = 0.50 

Fig.5 Caustic and critical lines in neighborhood of |Q| = 1/24 at which beak-to-beak 
appears on F-axis. 

Fig.5a \Q\ = 0.03 
Fig.5b \Q\ = 1/24 ~ 0.04 
Fig.5c \Q\ = 0.05 

Fig.6 Caustic and critical lines in neighborhood of |Q| = Qi = 1/8 at which butterfly 
appears. Point P in Fig.6b is the point where butterfly appears. 

Fig.6a \Q\ = 0.12 
Fig.6b \Q\ = 1/8 
Fig.6c \Q\ = 0.13 

Fig.7 Caustic and critical lines in neighborhood of |Q| = Q 2 = which 

swallowtail appears. 


Fig.Ta \Q\ = 0.17 
Fig.Tb \Q\ = g2 ~ 0.1745 
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Fig.Tc \Q\ = 0.18 


Fig.8 Shapes of the images, critical and caustic line for the multiple quasar PG1115+080 
from our calculation. The lensing galaxy is at coordinate origin. Observed image 
positions are denoted by crosses. 

Table 1 Observed image positions and their amplifications of PG1115+080. 
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image 

a;i(sec) 

a;2(sec) 

/i//iC 


-T'.27 

-2". 08 

3.22 

^2 

-1".44 

-T'.62 

2.49 

B 

0".39 

-T'.95 

0.64 

C 

o 

o 

o 

O 

o 

o 

1.00 

lens 

-0".33 

-T'.35 



Table 1: Observational data 
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